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Abstract 

We study numerically self-trapped (polaron) states of quasiparticles (electrons, holes or 
excitons) in a deformable nanotube formed by a hexagonal lattice, wrapped into a cylin- 
der (carbon- and boron nitride-type nanotube structures). We present a Hamiltonian for 
such a system taking into account an electron-phonon interaction, and determine conditions 
under which the lowest energy states are polarons. We compute a large class of numerical 
solutions of this model for a wide range of the parameters. We show that at not too strong 
electron-phonon coupling, the system admits ring-like localized solutions wrapped around the 
nanotube (the charge carrier is localized along the nanotube axis and uniformly distributed 
with respect to the azimuthal coordinate). At stronger coupling, solutions are localized on 
very few lattice sites in both directions of the nanotube. The transition from one type solu- 
tion to the other one depends on the diameter of the nanotube. We show that for the values 
of the carbon nanotube parameters, the polarons have a ring-like structure wrapped around 
the nanotube with a profile resembling that of the nonlinear Schrodinger soliton. 

1 Introduction 

Low-dimensional systems have been studied intensively during the last decades. Some systems 
have been found to possess unusual properties which turn out to be important for their prac- 
tical applications in microelectronics and nanotechnologies. Conducting polymers, quasi-one- 
dimensional organic and inorganic compounds are good examples of such systems [1]. Recently, 
nanotubes have become the subject of many intensive experimental and theoretical studies. In 
particular, carbon [2, 3, 4] and boron nitride [5] nanotubes have been shown to be useful for 
miniaturised electronic, mechanical, electromechanical and optoelectronic devices. Some of the 
intriguing physical properties of such systems are due to the electron-phonon interactions which 
are particularly important in low-dimensional systems. 

Although low-dimensional systems have been studied mainly within simplified single-chain 
(single-band) models resulting in the theoretical predictions of many striking phenomena [6, 7, 8], 
many of them, such as Peierls transitions, formation of charge density waves, soliton states, etc 
have been verified experimentally in real substances. However, taking into account the more 
realistic complex structure of these substances is quite important for the proper description of 
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their properties and can result in qualitatively new features. This is particularly true in the 
case of nanotubes who possess series of energy bands which are determined by one-dimensional 
(ID) energy dispersion relations on the wave vector k along the axis of the nanotube. It is well 
known that sufficiently long single wall carbon nanotubes (SWNT), in which half of the energy 
band states are occupied, possess properties of ID-metals or semiconductors depending on their 
diameter and chirality [2, 3, 9, 10]. In fact, their properties can be modified in a controllable 
way by varying either the nanotube diameter and chirality or by doping nanotubes with some 
impurity atoms, molecules or compounds [11]. 

In ID systems the electron-phonon interactions can lead to the formation of stable self- 
trapped states with spontaneous symmetry-breaking distortions [8, 12]. Such soliton states are 
realized in quasi-lD organic and inorganic materials and can successfully describe various prop- 
erties of conjugated polymers in the framework of the Su-Schrieffer-Heeger model [12]. The 
possibility for the formation of states with a spontaneously broken symmetry in carbon nan- 
otubes has been discussed in [10, 13, 14]. In particular, large polarons (or solitons) in nanotubes 
have recently been considered in [15, 16]. In these models the long-wave approximation has been 
used for the states close to the Fermi-level, which is equivalent to the continuum approximation. 
In this way, however, the essential details of the crystallographic structure of the system are not 
taken into account. 

Here we study numerically the formation and properties of polaron states in zigzag nanotubes 
using the semi-empirical tight-binding model in which we consider only the nearest-neighbour 
hopping interactions [2]. The merit of this method is its simplicity, in spite of which it often 
gives a quantitative agreement with the experimental results. This has been demonstrated in [12] 
for conjugated polymers and for carbon nanotubes in [2, 3, 9, 10]. Moreover, we consider only 
the ground state of a single electron in a hexagonal nanotube. Such systems are already quite 
complicated and though such states in ID chains have been studied intensively both analytically 
and numerically, little is known about them in the case of nanotubes. 

2 The Model 

We consider a model in which a two-dimensional sheet of hexagonally bound equivalent atoms 
(graphene-like sheet) is rolled-up into a zigzag nanotube. Calling d the length of each side of 
the hexagons, R the radius of the nanotube, n the number of hexagonal cells wrapped around 
the nanotube, we define (cfr. Fig 1) 

a = 2vr/n, a = dy/S, b = d/2. (1) 



The site labelling that we use is shown in Fig.l and coincides with the conventional one 
[2]. Here j is an axial index, i = 1, ...,n is the index labelling the sites around the nanotube, 
and p = 1, 2 labels the two types of atoms (usually called A and B): those who have a nearest 
neighbour one site down (p = 0) and those who have a nearest neighbour one site 'up' (p = 1) 
- see Fig.l. Then the position of any nanotube lattice site, at its equilibrium, can be described 
by: 

Rl hp = Rsm((i + J -^)a)e x + i?cos((* + ] -^)a)e y + ^^de z . (2) 

In Fig. 2 we present a drawing of the horizontal cross-section of the nanotube. Note that 
the parameters R, a and a are connected by 

a = 4i?sin(^). (3) 

Next we define 

Uij, p = Uij :P + Sij iP + Vij, p (4) 
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Figure 1: Hexagonal lattice, j is assumed even. 
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as the displacement vector of the atoms from their equilibrium positions where Uij tP is tangent 
to the surface of the undeformed nanotube and perpendicular to the nanotube 
tangent to this surface and parallel to the nanotube axis, and Sij jP is normal to the surface of 
the nanotube. Then, using Cartesian coordinates, we have 

Ui,j, p = «ij, p (cos((i + (j + p)/2)a)e x - sin((i + (j + p)/2)a)e y ), 
Si,j, P = Sjj, P (sin((i + (j + p)/2)a)e x + cos((z + (j + p)/2)a)e y ), 

v i,j,p = v i,j,p e zj (5) 

where Uij tP , Si t j tP and Vij tP are the corresponding amplitudes of displacements. So the positions 
of the lattice sites are given by 

Ri,j, P = Ri,j :P + U i:j ,p. (6) 

We also define the following lattice vectors connecting the atom p) with its three nearest 
neighbours 5(i,j,p) with 8 = r,l,d for right (r), left (I) and down or up (d) neighbours: 

D8i,j, p = Rs(i,j, P ) — Ri,j,p- (7) 

When Ui,j tP = we add the upper index to all quantities to indicate their values at the 
equilibrium position. Hence D5ij tP = DS^j p + dSij tP with dSi t j tP = Us(i,j, P ) ~ ^i,j,p- Note that 

\Dr-jJ = \Dl° hJ J = \Ddlj = d. (8) 

The potential energy of the lattice distortion includes central forces, which depend only on 
the distance between two sites. In the case of small displacement, i.e when \d8i : j tP \ <C d, the 
distance between lattice sites is approximately given by: 

\DSij >p \^d + WS iJ>p , (9) 

where 

d8i j o ■ D5a a n , s 

wSijj = w d w (10) 

are the changes of distances between nearest neighbours due to the displacements of the sites. 

As is well known, the central forces between neighbouring sites do not provide lattice stability. 
Therefore, we introduce additional terms. First, assuming that the force between two atoms is 
not purely central, in addition to W8ij jP which are invariant under translations, we introduce 
also the quantities Q5ij P which describe relative shifts of neighbouring atoms in a graphene sheet 
and which will be needed to remove the non-physical zero modes of the phonon Hamiltonian. 
The explicit expressions for these quantities are: 

nr i,jfi = ^( cos (|)Kj,i - u i,j,o) + sin(^)(sij,i + sijflf) - y(fij,i - Uij,o), 

n kjfi = ^( c °s(^)(^,j,o - iM-ij,i) + sin(|)(s i _i, ji i + Sijflf) - ^(ui-ij,! - v iJfi ), 
Qdijfi = —Uij-i^i + Uijfi 

^ r i,j,l = ^ r i,j,0, ^k,j,l = ^k+l,j,0, ®*di,j,l = ^<k,j+l,0- (11) 



Next, as in [17, 18], we introduce also the force of bond-bending between the atoms induced 
by the deviation of the valence angles from their equilibrium values. For this we compute the 
solid angle spanned by the 3 lattice vectors located at a given site. Clearly, this angle is given 
by 

q _ (Dlj,j,p x D r i,j,p)-Ddj,j, p , , 

\ JJr i,j,p\\- L "'i,j,p\ \ J -JO'i,j,p\ 
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For small displacements we have 



Si,j,p ~ Sfj if) + ^Cij tP (13) 



where Sfj p = |sin(|) and 



\/3 . ,ck 

Ci,i,o = — sin(-)(2u iJi o - v iJtl - Vi- 1Jtl ) 
-cos(-)s ij _i i i + 3cos {-^jSijfl 

,3 , 0L n 5 0,0;,.. . 
+ (g cos ( 4 ) - 2 COS (4 + s m.i) 

+ sin (^)(2 cos (4) ~ ^Kj,! 

Cjj,i = — sin(-)(uy i0 + Ui+ij.0 - 2v iJ):l ) 

- cos(-)a i j + i,o + 3 cos (-)sij, 1 

.3 , CK , 5 o . . . . , 

+(g cos ( 4 ) - 2 COS ( 4 ))( Si+1 J'° + s m>°) 

+ sm(^)(^cos 2 (^) - l)(u mj - - «ij,o)- (14) 

Taking into account the harmonic potential terms responsible for the central i.e Vw, non- 
central, Vq, and the bond-bending, Vq forces, and using the approximation of the nearest- 
neighbour interaction, we define the phonon Hamiltonian as: 

it _ 1 \ " f Pi,j,p 1 %hP 1 T ij,P 1 j, r< 2 
H P h - 2MIW" + ~W + ~M + kcU ^ 
h3,p 

+ k [Wrl hp + Wll hp + Wdl hp + Ori H , + ^ p + tl^J) . (15) 

Here Pij tP , Qi,j,p and fi,j,p are the momenta, canonically conjugate, respectively, to the displace- 
ments Uij jP , Sij tP and Vij tP . M is the atom mass, k is the elasticity constant that characterises 
the central force, and k c is a characteristic constant of the bond-bending force. 

Next we assume that each atom in the hexagonal lattice has a well isolated non-degenerate 
electron level £q (-zr-orbital in the case of a carbon nanotube) and we take into account the 
hopping interaction, J, of each atom only with its three nearest neighbours. 

The electron-phonon interaction has different sources [10, 19, 13, 17, 20]. Often, the phonon 
modulation of the hopping interaction is invoked which, in the linear approximation with respect 
to the displacements, gives 

J(i,j,p);8(i,j,p) = J - G2W5ij tP , (16) 

thus leading to the electron phonon interaction proportional to GV 

In general, as has been shown in [17], in addition to the conventional modulation of the 
hopping interaction one should also take into account the interaction between electrons and ions 
through the Coulomb potential which is modulated by the lattice vibrations. 

The fields of neighbouring atoms then alter the local electron energy, £ij tP , and so, in the 
same linear approximation, we can write 

£i,j, P = £o + Xi (Wrij, p + WUj,p + Wdij, p ) + X2 C iJiP , (17) 

where \i characterises the influence of the distance and X2 the influence of the bond-bending 
deviations. 
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The constants G2, Xi an d X2 determine the strength of the electron-phonon coupling. Their 
calculation would be the subject of microscopic considerations. The estimate of G2 for carbon 
nanotubes was obtained in [21]. The coupling, through the Coulomb potential, provides an 
input into all these constants. Moreover, the curvature of the lattice is an important factor: 
the curvature-induced charge polarisation, present in carbon nanotubes, results in the vibration 
modulated dipole component of the Coulomb potential. As it was indicated in [22], where this 
fact had been taken into account through the Lydann-Sachs- Teller relation, such an electron- 
phonon interaction is stronger that the one due to the deformational potential. In this paper we 
treat the constants G2, Xi an d X2 as parameters and we investigate the ground self-trapped state 
of an electron depending on the value of these parameters. The strength of the electron-phonon 
interaction plays a key role in many effects. When the coupling constant of this interaction is 
strong enough it can lead to the self-trapping of a quasiparticle. By 'quasiparticle' we denote 
an electron, a hole or an exciton in a nanotube. 

The self-trapped states of quasiparticles are usually described in the adiabatic approximation 
which is equivalent to the semiclassical consideration in which the vibrational subsystem is 
described as a classical one. Therefore, the state of a quasiparticle in such a model is described 
by the Hamiltonian functional 



h3,P 



+ Xl l^ij.pl 2 ( Wr i,j,P + W h,j,p + Wd i,j,p) + X2 IVijJ 2 Ci,j,a 

+ 9 G 2[« j , p ^r(i,j,p) + ti(iJ,p)<Pij,p)Wr iijtP + (<Pij,p<Pl{ij,p) + tf(iJ,p)<Pij,p)Wlij,p 



+ 



(18) 



Here H p h is given by (15) and fij, p is the amplitude of the probability of quasiparticle being 
located on the site (i,j,p). 

From (18), after expanding the W's, fi's and C's, we can derive the following static equations 
for the functions <fi,j, P , u, v and s: 







(W + £o)^P ~ J (<Pr + <Pl + <Pd) +Xl<P (W r + Wi + W d )+X2<pC 

+ G 2 [f r W r + wWi + <p d W d ], 



(19) 



1— ( X IX 

>/3cos(-)(Wj - W r ) + cos(-)(fy - n r ) + 2Q d 



sin(|)(^cos 2 (|)-l)(Q-a) 



+ — C0S(-)G 2 (V3 <Pl + <Pi<p- V <Pr + <Pr t P), 



,a 5 
4^2 



•4- 



(20) 



(2w d - w r - Wi) + V3(n r + n t 



^sin(|)(2C + C r + Q) 



+ xi 2 ( 2 l^d| 2 " WA 2 ~ Wl\ 2 ) + X2(2|^| 2 + Wr? + Wl\ 2 ){— sin(")) 



+ ^G2(-^>*fr ~ f*f ~ ~ <P*<P + ty*<Pd + 2( P d <P), 



(21) 



V3sm(^)(Wr + Wl) + sin(^)(ft r + ft,) 
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+ k. 



(| cos (l) " \ «* 3 (f + C/) - cos(|)C d + 3 cos 3 (^)C 



+ Xl ^ s in(|)(2|^| 2 + |^ r | 2 + |^| 2 ) 



+ X2((g cos(^) - - cos 3 (|))(|v9 r | 2 + |^| 2 ) - co S (|)|^| 2 + 3cos 3 (^)M 2 ) 



+ G 2 —sm(-)((p*ip r + (p* r ip + ip*(pi + tp* l (p). 



(22) 



Here we have used the notation <p, ip r ,ipi, ip^ for the functions on the lattice site {i, j, p} and on 
its corresponding 3 neighbours. 



3 Numerical results 

In this section we describe various numerical stationary solutions of our equations. Note that 
our equations possess several parameters: J, k, k c , G 2 , Xi an d the lattice spacing d. For the 
numerical calculations we need dimensionless quantities and so we perform the following change 
of variables: 

u = IU; v = IV; s = IS; k = -kJ; k c = -prJ; 

Xi = jtJ; X2 = jj-J; G 2 = -^J; £ = £ J. (23) 

First we choose to measure the energy in units of J and so we put J = 1. Then So is dimen- 
sionless. Next we choose the unit of length I to be such that k = 1. From the expressions above 
we see that this corresponds to / = \J J/k. For the carbon nanotube the parameters take the 
following values J « 2AeV [2, 20, 4], d = 1.42 x l(T 10 m [2], k = 36.5 x W 4 dyn/cm [2]. Using 
these values, we get I = (3.24 - 3.7) x KT n m. Thus, the length scale we use is smaller than 
the lattice length d « 4Z. For carbon nanotubes the dimensionless electron-phonon coupling 
constant G = G 2 l/J ~ 0.6 - 0.74. 

We have derived numerical solutions of our model for all values of n between 3 and 20 for 
several ranges of other parameters. In all our calculations we have taken k = 1, k c = 0.2. Then 
fixing the values of two of the three parameters X±, X 2 and G we have varied the value of the 
3rd one between and 5, and xi and X2 mostly between and 1. We have also set the lattice 
parameter a = 1 and, in our discussions, we will refer to the main axis of the nanotube as y 
while x will run along the circumferences of the nanotube. 

Of course, we have used periodic boundary conditions along the x axis, while for the ends 
of the nanotube we have used free boundary conditions. This allows not only for a global 
translation of the nanotube but also for a dilation or contraction of the lattice along its axis. To 
exclude global translation of the lattice, we have pinned a single lattice site, at the bottom of 
the lattice, using a quadratic potential. Using such a potential, rather than fixing the position 
to a fixed value, is smoother from the numerical point of view. 



3.1 Numerical solutions 

We have found that, depending on the values of the parameters, the system of equations (19-22) 
possesses several types of solutions. Some of them, sometimes, coexist at the same values of the 
coupling constant. Some of these solutions are shown in Figs. 3-9 where we have plotted the 
hexagonal lattice unfolded onto a plane. The periodic boundary conditions imply that the left 
and right hand sides of the lattice are connected together (the lines joining these sites are not 
drawn in our figures, they are implicit). The position of each site is given by its equilibrium 
value translated by the deformations u and v at that lattice site. The s field, on the other hand, 
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is represented by an arrow, on the same scale, which points downwards if the s field points 
towards the centre of the tube, i.e. if the tube is pinched, or upwards if the field points away 
from the tube, i.e. if the tube is dilated. 

Most of the solutions we present here correspond to the n = 8 case but the solutions for 
other values of n look very much the same. Occasionally we will add some comments on this 
dependence. In the figures we only present a small section of the nanotube we have taken; for 
clarity of the picture we do show the full width of the tube, along its periodic direction, but 
we restrict ourselves to the section of the tube along its axis where the soliton is located. The 
solutions were actually computed on a grid ranging from —50 to +50. 

In our plots we have used the dimensionless units described above, which correspond to 
displacements 4 times larger than the real displacements of the nanotube, except for Figure 3. a 
where we had to amplify the dimensionless units by a factor of 10 and Figure 4.b where we have 
used the physical units because the displacements in this case were already sufficiently large. In 
the figures, each lattice site is represented by a circle, the inside of which, is painted in shades 
of black describing the value of the polaron probability density |c/?| 2 . A white inside of the circle 
corresponds to a very small density, while black is used for the largest density. The actual shade 
scale is given on the right hand side of each figure. Note that the axis of the nanotube runs 
vertically in our figures. 




Figure 3: Solution for n = 8, X\ = 0.6, X2 = 0.2 and a) G = 1. b) G = 1.5. The deformed 
lattice is plotted with small circles at the vertices representing the electron density (see the grey 
scale on the right; white corresponding to \tp\ 2 = 0). The arrows at each vertex represent the 
s field; a down/up arrow represents a displacement towards/away from the centre of the tube. 
The displacements are plotted in dimensionless units for G = 1.5 and amplified by a factor of 
10 for G = 1. 

When the electron-phonon coupling constants are small the solutions are invariably ID bell- 
shape solitons wrapped around the nanotube. This is seen in Fig. 3. a where the amplitudes 
of the displacements have been amplified 40 times compared to their physical values. In what 
follows, we name these solutions - type I solutions. We see that for type I solutions the polaron 
probability distribution looks like a ID Davydov soliton along the axis of the nanotube which 
is spread uniformly around it. For our choice of parameters the nanotube is contracted at the 
position of the soliton and the tube itself is also pinched. Note that if we choose the 3 electron- 
phonon coupling constants with the opposite signs, the displacements fields u, v and s will 
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change signs and the nanotube will become extended and inflated. 

Note also that the soliton drags the whole nanotube towards itself. In other words, the 
contraction of the lattice at the soliton position is not compensated by an appropriate dilation 
elsewhere. 

When we increase the electron-phonon coupling constants above some critical value, the 
radially symmetrical solutions do not correspond to the minimum value of energy and the polaron 
localises on a few sites. We show some typical solutions in Figs. 3.b, 4 and 9. The solution at 
X\ = 0.6, X2 = 0.2, G = 1.5 is shown in Fig. 3.b): the polaron probability density reaches a 
maximum on a single lattice site, but the polaron is spread over several lattice sites around the 
position of its maximum density. The lattice is slightly contracted but the lattice sites, where 
the polaron is localised, are slightly distorted into the nanotube. We name such solutions - 'the 
type II solutions'. 




Figure 4: Solution for n = 8,X 1 = 0.6, X 2 = 0.2 and a) G = 2 b) G = 5. The deformed lattice 
is plotted with small circles at the vertices representing the electron density (see the grey scale 
on the right; white corresponding to |i^| 2 = 0). The arrows at each vertex represent the s field; 
a down/up arrow represents a displacement towards/away from the centre of the tube. The 
displacements are plotted in dimensionless units for G = 2 and in physical units (dimensionless 
units divided by 4) for G = 5. 

In Fig. 4. a, we present a solution which is also localised on a few sites but for which the 
maximum of the polaron density is evenly distributed between 2 lattice sites. This solution 
corresponds to the same values of parameters as figure 3.b except that G = 2. In this case, the 
lattice is also contracted, the polaron is more localised, and the lattice sites where the polaron 
is localised are also moved into the nanotube. Notice that the two sites at which the polaron 
density reaches its maximum have moved to be very close to each other. We name such solutions 
- 'the type III solutions'. 

In Fig. 4.c, we present a solution corresponding to the same values of the parameters as 
figures 3.b and 4. a but for which G = 5. This time the polaron is localised mainly on 3 lattice 
sites: it has a probability density of about 0.5 on one site and of 0.25 on either of the two 
neighbouring sites. The lattice deformation in this case is very large and so we had to scale it 
down to the physical scale. Once again, the lattice is contracted at the position of the polaron 
and the lattice sites where the polaron is located have all moved into the nanotube. We name 
these solutions - 'the type IV solutions'. 
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The four classes of solutions we have presented so far, correspond to the same family of 
solutions obtained when one varies G and takes X\ = 0.6 and X2 = 0.2. In Figs. 5-7 we present 
the maxima of the polaron probability densities (D = M.axuj iP \ \<Pij,p\ 2 ) and the eigenvalues A 
of these solutions as a function of G. We present the data for the cases of n = 20, 8 and 3. 
The solutions shown in figures 3 and 4 correspond to the 4 types of configurations presented in 
figures 6. In particular we note that the transitions between the different types of solutions are 
always sharp and that sometimes there is an overlap: two different types of solutions coexist for 
the same values of the parameters. 

We also note that the diameter of the nanotube does not affect very much the types of 
solutions that exist, except for the location of the transition between the different types of 
solutions. On the other hand, the binding energies of the solutions are larger for smaller values 
of n. 



D 




Figure 5: The maximum of probability, D = Max^- p) \tyi,j,p\ 2 and the eigenvalue A for X\ = 0.6, 
X2 = 0.2 as a function of G for n = 20. 

In Fig. 8 we present the values of G for which the various types of solutions exist as a 
function of n. We see that there is an overlap between solutions of type I and II and that this 
overlap is larger for larger values of n. 

In Fig. 9. a, X\ = 2, X2 = 0, G = 0.6 we present a solution which appears to be rather 
different from the previous ones. This time, the polaron is localised strongly on 4 sites, but with 
a density of about 0.55 at the central site. Note that the deformation of the lattice is nearly 
symmetric this time. A slight asymmetry is induced by the bending of the nanotube. We name 
these solutions - 'the type V solutions'. 

Finally in Fig. 9.b we present a solution that is strongly localised on a single site. Note that 
this time the lattice is very strongly deformed. We name such solutions - 'the type VI solutions'. 

3.2 Solutions for X\ = 0.6, G = 0.6 and varying X 2 

In Fig. 10 we present the plots of the values of the maximum density and of the eigenvalues of 
the solutions for X\ = 0.6, G = 0.6 and n = 20 as a function of X%. We note that in this case 
all solutions fall into two classes seen before. The solutions with a small value of the maximum 
density are of type I while the highly peaked ones are of type VI. The solutions look very much 
the same for other values of n except that the ranges of parameters for which the solutions exist 
do depend on n. 
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Figure 6: The maximum of probability, D = Max^- p) \tpij tP \ 2 and the eigenvalue A for X\ = 0.6, 
X 2 = 0.2 clS Si function of G for n = 8. 




Figure 7: The maximum of probability, D = Max^j^l^l 2 and the eigenvalue A for X\ = 0.6, 
X2 = 0.2 as a function of G for n = 3. 
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Figure 8: Domain of solutions as a function of n for solutions with X\ = 0.6, X 2 = 0.2. o: upper 
value of G for ring like solution (type I), x: lower value of G for solutions of type II. • : Value 
of G at the transition between solutions of type II and type III. + : Value of G at the transition 
between solutions of type III and type IV. 




Figure 9: Solution for n = 8, G = 0.6 and a) XI = 2, X 2 = b) XI = 0.6, X 2 = 1.5. 
The deformed lattice is plotted with small circles at the vertices representing the electron density 
(see the grey scale on the right; white corresponding to \(p\ 2 = 0). The arrows at each vertex 
represent the s field; a down/up arrow represents a displacement towards/away from the centre 
of the tube. The displacements are plotted in dimensionless units. 
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Figure 10: The maximum of probability, D = Max( ij)/9 ) |9?i.j,p| 2 and the eigenvalue A for X\ = 
0.6, G = 0.6 as a function of Xi for n = 20. 

In Fig. 11 we present the regions of parameters for which the solutions of type I and type 
VI exist as a function of n and Xi- Solutions of type I exist for all values of X2 below the circles 
while those of type VI exist for all values of X2 above the crosses. 

We see that there is a large overlap of the regions in which the solutions of the type I and 
VI coexist and that the region of this overlap grows with n. 

3.3 Solutions for X 2 = 0, G = 0.6 and various values of X\ 

In Fig. 12 we present the plots of the maximum density and of the eigenvalue of the solutions 
for X2 = 0, G = 0.6 and n = 20 as a function of X\. Again we note that there are two types 
of solutions in this case. The solutions with a small density maximum are of type I while the 
highly peaked ones are of type V. This time, however, there is very little overlap between these 
two classes of solutions. 

The solutions look very much the same for other values of n, except that, again, the range 
of parameters where the concrete solutions exist does depend on n. For a given value of the 
parameters, the solutions are also more bound when n is small. 

From this numerical study and what we know about the physical values of the parameters 
in our model we can conclude that the type I solutions exist in all real nanotubes. If X\ and 
X2 are large, then strongly localised solutions of other types may exist as well. 

4 Conclusions 

In this paper we have determined the conditions for the existence of polaron (narrow soliton) 
states in hexagonal nanotube systems. They are given by the solutions of the equations minimis- 
ing the total Hamiltonian of the system. We have also studied the properties of these solutions 
showing that in some cases the system possesses two or more different states very close in en- 
ergy. Some of these states are localised on few lattice sites, some are more spread out. We 
have determined the boundaries between various states and have found that in some cases these 
boundaries are very sharp. The corresponding critical values of the coupling constants depend 
on the diameter of a nanotube. We have also shown that the solutions corresponding to the 
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Figure 11: Domain of solutions as a function of n for solutions with X\ = 0.6, G = 0.6. o: upper 
value of X2 for ring-like solutions (type I), x: lower value of X2 for solutions of type VI. 




Figure 12: The maximum of probability, D = Max^jpj l^ijjp) 2 an d the eigenvalue for X2 = 0, 
G = 0.6 as a function of X\ for n = 20. 



14 



physical parameter values of the carbon nanotube are ring-like solutions wrapped around the 
nanotube with a profile resembling that of the nonlinear Schrodinger soliton, i.e. similar to a 
ID soliton (comp. [15]). 

Our numerical results do not apply straightforwardly to real carbon nanotubes with one 
electron per carbon atom when half of the electronic band states are occupied. Nevertheless, 
our results may be applicable to doped nanotubes at doping levels such that the concentration 
of charge carriers is sufficiently low so that the distances between self-trapped carriers exceed 
their localisation width. As has been shown in [23], in the case of the self-trapping, the Fermi 
statistics manifests itself by the appearance of the spatial repulsion of bipolarons. Two electrons 
with opposite spins can occupy the same place forming a bipolaron state. Two electrons with 
parallel spins or four electrons, in general, can create two self-trapped states. These states 
are characterised by a Fermi repulsion between them which decreases exponentially with the 
increases of the distance between them [23]. When the average distance between (bi)polarons 
is larger than the localisation region, the interaction between them can be neglected, and the 
single-(bi)polaron model, used in the present paper, becomes justified. 
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